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1. Introduction and motivation 

The theoretical investigation of gauge theories in the \nrge.-N limit has been a very active 
research field for almost forty years: the most relevant early works, including the seminal paper by 
't Hooft [1], are collected in ref. [2]; in addition, several lecture notes and reviews are also available 
in the literature, see, e.g., refs. [3, 4]. A new review article — covering, in paiticulai^, most of the 
lattice studies of large-A^ gauge theories — ^has recently appeared [5]; the continued interest of the 
lattice community in this topic is confirmed by the large number of contributions in the parallel 
sessions of this conference [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. 

1.1 QCD in the 't Hooft limit 

Large-A^ gauge theories are generalizations of QCD in which the number of color charges A'^ 
is taken to be infinite. In order for this limit to make sense, at least perturbatively, at the same time 
the coupling g has to be taken to zero, holding the 't Hooft coupling A = g^N fixed. Keeping also 
the number of quark flavors Uf fixed yields the so-called 't Hooft limit [1].^ 

This double limit is characterized by interesting properties, which dramatically simplify the 
theory. In particular, the perturbative dynamics is dominated by a special class of Feynman di- 
agrams: counting the powers of g and keeping track of the number of independent fundamental 
color indices appearing in a Feynman diagram, through the double-line notation — which repre- 
sents propagators of quarks (gluons) in the fundamental (adjoint) representation of the gauge group 
by single (double) lines — , it is easy to see that the contributions with the largest power of A'^ come 
from planar diagrams {i.e., diagrams which, in this double-line notation, can be drawn on a plane 
without any crossing lines) without dynamical quark loops: see fig. 1.1 for an example. 

More generally, ampUtudes for physical processes can be expressed in terms of double series, 
not only in powers of A, but also in powers of 1/A'^. The latter expansion has a topological nature, 
as it con^esponds to an expansion in classes of diagrams which can be drawn on Riemann surfaces 
of different topology, i.e. with a different number of "handles" {h) and boundaries {b): 

oo oo 

^= £ A^'-'"-''Ic(;,,),„A". (1.1) 

h,b=Q n=Q 

Interestingly, a similar topological expansion is also found in string theory, upon replacing \/N 
with the string coupling gg. This observation led to speculations that string theory could actually be 
a reformulation of \wg&-N gauge theory already during the 1970's. From a more modern perspec- 
tive, this correspondence is also expected to hold in the conjectured duality relating gauge theories 
and string theory [19, 20, 21] : in the large-A'^ limit, loop effects on the string side become negligible 
(see also the plenary contribution by Hanada [22]). 

1.2 A wealtii of plienomenological implications from large-A'^ counting rules 

The large-A'^ counting rules introduced above have many phenomenologically interesting im- 
plications. In particular, if one assumes that QCD in the 't Hooft limit is a confining theory, it is 

'a different type of large-W limit (Veneziano limit) is obtained, if one also takes the number of quark flavors 
to infinity, holding iif/N fixed [18]. As this limit is generally less easy to study than the 't Hooft limit, the present 
discussion will be mostly focused on the latter. 
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Figure 1: Top: Different types of Feynman diagrams contributing to the gluon propagator at three-loop 
order: in the 't Hooft Hmit, the number of independent fundamental color indices is largest for the diagram 
on the left panel, which only includes planar gluon loops. Replacing a gluon loop with a quark loop (central 
panel) reduces the power of by one. Finally, the diagram on the right panel (including only gluon loops, but 
with a non-planar topology) is (9(1). As a consequence, in the 't Hooft limit the diagrams on the central and 
right panels are subleading with respect to the one on the left panel. Bottom: Two diagrams contributing to 
the propagation of a meson (denoted by a blob): the diagram on the right, in which the process goes through 
an intermediate stage involving only two virtual gluons, is suppressed by a factor of 1 /N with respect to the 
one on the left, as expected according to the OZI rule. 

possible to deduce that the spectrum consists of an infinite number of infinitely narrow mesons and 
glueballs, with masses 0(1). Their interactions are suppressed by powers of 1/^/N, so that large-A^ 
QCD is a theory of stable, weakly interacting hadrons. Secondly, exotic states (e.g., tetraquarks, 
molecules, et c.) are absent. Thirdly, the OZI rule [23, 24, 25] is exact — see fig. 1. 1 for an example. 

Another implication of the large-A^ counting rules is that loop effects in the effective chiral 
Lagrangian ai^e suppressed, so that it can be studied in the tree-level approximation. Moreover, the 
axial anomaly turns out to be proportional to 1 /N, and hence suppressed in 't Hooft's limit. 

Baryons can be interpreted as the solitons of large-A'^ QCD, with masses 0{N) [26]; further- 
more, by imposing certain consistency conditions related to unitarity (see refs. [27, 28] and refer- 
ences therein), it is possible to derive a systematic expansion in powers of 1 /N for quantities such 
as baryon-meson couplings, baryon masses, magnetic moments, et c. 

The large-A'^ limit might also have implications for the QCD phase diagram: in particular, 
McLerran and Pisarski conjectured the existence of a new, exotic state of matter {quarkyonic mat- 
ter) at large densities [29] — see also ref. [30] for a related lattice model. 

Finally, the large-A'^ limit also has a number of interesting implications relevant in the high- 
energy domain (evolution equations, hadronic cross-sections, paiton distributions and structui'e 
functions, large-A'^ Standard Model, et c.) [31], that will not be discussed further here. 

1.3 Factorization, volume reduction and large-A'^ equivalences 

Besides interesting phenomenological aspects, \wge.-N gauge theories are also characterized 
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by many intriguing properties at the fundamental level. In particular, the large-A'^ counting rules im- 
ply that vacuum expectation values (vev's) of products of gauge-invariant operators are dominated 
by disconnected contributions. This immediately leads to the factorization of vev's of physical 
operators, up to 0{\/N) corrections: 

(^1^2> = (^l)(^2>+0(l/^), (1.2) 

which reveals an analogy with the classical limit of a quantum theory (with 1/A^ playing the role of 
fi). In fact, the analogy can be made explicit, by constructing an appropriate set of coherent states, 
and a "classical" Hamiltonian [32]. 

Factorization has an interesting consequence: the Schwinger-Dyson equations satisfied by 
Wilson loops in the large-A'^ theory on the lattice are independent of the physical hypervolume of 
the system, provided center symmetry is unbroken [33]: this is the so-called Eguchi-Kawai (EK) 
volume reduction. In principle, this property would allow one to study the large-A'^ theory in ar- 
bitrarily small volumes, either by analytical techniques (reducing the original theory to a matrix 
model), or by numerical simulations on a single-site lattice. However, it is well-known that center 
symmetry does get broken in a small volume in the continuum limit: this can already be seen at 
the perturbative level, for all D > 2. In order to preserve center symmetry, various fixes have been 
proposed, since the 1980's: for example, in the quenched EK model [34], one studies the dynam- 
ics of the single-site model for a fixed set of eigenvalues of the link variables along the various 
directions, and then averages over a center-symmetric distribution for the eigenvalues. However, 
this method has recently been shown to fail [35], due to the fact that the quenching prescription 
fixes the eigenvalues of the link matrices in the four directions only up to cyclic permutations, 
and dynamical fluctuations lead to non-trivial correlations among the eigenvalues along different 
directions. Another approach to preserve center symmetry in the reduced EK model is based on 
imposing twisted boundary conditions [36, 37]: interestingly, this approach can also be used for 
a non-perturbative definition of field theories in non-commutative spaces [38, 39, 40, 41, 42] (for 
an alternative regularization of these theories, see, e.g., refs. [43, 44] and references therein). Vol- 
ume independence in the twisted EK model holds both at strong coupling and in the perturbative 
regime, but, at least for the simplest definition of the twist, it has been found to fail at intermediate 
couplings, in a range which appears to increase when grows [45, 46]. However, a couple of years 
ago, the authors who originally proposed the twisted EK model suggested a new formulation of the 
twist [47], which they are currently studying numerically [48]. As an example, the plot on the left 
panel of fig. 1.3, taken from ref. [48], shows how the extrapolation of results for the string tension 
obtained from simulations in large volume (at moderately large A'^) compares with the result from 
a single-site simulation in the new version of the twisted model, at a much larger value of A'^: the 
results appeal- to be in nice agreement, and motivate further studies of this model. 

Another possibility to preserve center symmetry in EK models is based on the inclusion of 
dynamical adjoint fermions (with periodic boundary conditions in all directions): this idea was 
initially proposed in ref. [49], and has since been studied both analytically and numerically by a 
number of authors [13, 16, 50, 51, 52, 53, 54, 55, 56, 57]. Recent results indicate that, indeed, 
EK volume reduction with adjoint Dirac fermions works as expected, both with «/ = 1 and = 2 
flavors: see, e.g., the plot on the right panel of fig. 1.3, taken from ref. [54], which shows the center 
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Figure 2: Left panel: A comparison of results for the string tension obtained from simulations in a large 
volume (red symbols) and their linear extrapolation in 1 /N'^ (blue line), with the result for the same quantity 
from a single-site simulation in the twisted EK model (green symbol), from ref. [48]. Right panel: Numerical 
simulations of large-A^ gauge theories with dynamical fermions in the adjoint representation of the gauge 
group in a small volume show evidence that center symmetry is unbroken in a funnel-shaped region that 
extends to rather large quark masses, and remains of finite width in the continuum limit. The figure, taken 
from ref. [54], is a sketch of the phase diagram based on the results of simulations with nf — 2 adjoint flavors. 

symmetry realizations in the rif = 2 model, as a function of the hopping parameter K and of the 
bare lattice coupling b = 

Another way to enforce center symmetry in EK models is based on double-trace deforma- 
tions [58]: one modifies the usual Yang-lVIills (YIVI) action adding (products of) traces of Polyakov 
loops, with positive coefficients. This suppresses center-symmetry breaking configurations, at the 
cost of 0{l/N^) corrections to the observables: 

Sym -^S,u + ^Y. L HL"{x))\^ . (1.3) 

Related ideas have also been discussed in the context of SU(A'^) YIVI theory at finite temperature [59, 
60]. Dedicated numerical algorithms to study the EK model with double-trace deformation have 
been devised [61], and preliminary investigations are currently under way [62]. A nice feature of 
this approach is that one can reduce only one (or a few) direction(s), while keeping the others large. 

Finally, Neuberger and collaborators proposed the partial reduction approach to EK [63, 64]: 
the idea is to simulate the large-A/^ theory in lattices which are small, but still larger than the critical 
size corresponding to the inverse of the deconfinement temperature T^. As an example of results 
obtained in this approach, in ref. [65] the confining potential was computed up to distances equal 
to 9 lattice spacings, from simulations on a lattice of linear size L = 6a. 

Volume reduction and volume independence in large-A'^ gauge theories can be interpreted as an 
example of "orbifold" equivalence [66, 67], namely as a correspondence based on projections under 
some discrete subgroup of the global symmetries of two different theories [49, 68, 69, 70, 71, 72]: 
under the assumption that the discrete symmetry used in this projection is not spontaneously bro- 
ken, the vev's and correlation functions of invariant (or "neutral") sectors of observables in the 
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Figure 3: Left panel: Torelon spectram in SU(5) YM theory, taken from ref. [91], in comparison with the 
predictions fi^om the Nambu-Goto elfective string model. Central panel; Large-A^ YM glueball spectrum, 
taken from ref. [95]; the figure also shows the results for the ground state and first excitation in the J^*- = 0++ 
channel, and for the 2++ ground state taken from ref. [96]. Right panel: Dependence of the mass anomalous 
dimension in QCD^v, with different numbers of colors and 2 flavors of dynamical fermions in the two-index 
symmetric representation of the gauge group, on the re normalized 't Hooft coupling in the SF scheme, from 
ref. [104]. The dashed green line denotes the leading-order perturbative prediction in the large-A^ limit. 

original ("parent") and projected ("daughter") theories are equal — up to a trivial rescaling of cou- 
pling constants and volume factors. Such orbifold equivalences do not relate only theories defined 
in different volumes, but also theories with different field content: for example, the orientifold 
planar equivalence [73, 74] (investigated on the lattice in ref. [75]) can be interpreted as a coitc- 
spondence between two different daughter theories obtained by orbifold projections from a com- 
mon parent theory [72]. Finally, orbifold projections are also relevant for lattice formulations of 
supersymmetry [76] (see also refs. [77, 78, 79, 80] for work on related topics). 

2. A selection of recent results 

In this section, we discuss a selection of recent lattice results obtained from simulations of 
\nrge.-N gauge theories (in a large volume). 

2.1 Results in four spacetime dimensions 

As mentioned in subsec. 1 .2, many phenomenologically interesting implications of the \wge.-N 
counting rules are derived under the assumption that QCD (or YIVI) is confining in the large-A'^ limit. 
Testing the correctness of this assumption at the non-perturbative level, via lattice simulations, 
was one of the first problems to be addressed, and by now we know that, indeed, SU(A'^) gauge 
theories are confining in the 't Hooft limit (see, e.g., refs. [81, 82]). Lattice studies also show that 
confining flux tubes can be modelled quite accurately as Nambu-Goto strings (a fact generically 
observed in confining theories [83, 84, 85, 86, 87, 88, 89, 90]). As an example, the plot on the 
left panel of fig. 3, taken from ref. [91], shows a comparison of the torelon spectrum in SU(5) YM 
theory with the predictions from the Nambu-Goto model. Other works studying the effective string 
picture at large A'^ include those by Lucini and Teper [92], by Lohmayer and Neuberger [93] and by 
IVlykkanen [94]. 
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Figure 4: Left panel: Extrapolation of the ratio of the critical deconfinement temperature over the square 
root of the string tension to the large-A^ limit, from ref. [116]. Right panel: The pressure (p) and energy 
density (e), normalized to their respective SB limits, as a function of T /Tc, in YM theories with = 3, 4 
and 6 colors, from ref. [11 4]. 



As for the YM hadron spectrum, glueball masses turn out to have a smooth dependence on A^. 
The central panel of fig. 3, taken from ref. [95], shows the results extrapolated to the 't Hooft limit 
for the masses of ground-state and excited glueballs, in several different channels. The results from 
an earlier study [96] are also shown. 

Another interesting issue deserving non-perturbative investigation is the dependence of the 
large-A^ theory on the 't Hooft coupling. As mentioned in subsec. 1.1, the expectation that A is the 
appropriate coupling for the 't Hooft limit of the theory is based on perturbative arguments. Hence, 
one may wonder, whether this is borne out non-perturbatively. Lattice simulations do confirm that 
this is the case. This can already be seen from the bare lattice coupling: ref. [97] studied how 
the tree-level improved bare lattice 't Hooft coupling runs with the lattice spacing a (determined 
non-perturbatively from computations of the string tension) in YM theories with a different number 
of colors: the collapse of data obtained from simulations at different values of A'^ showed clearly 
that A is an appropriate coupling for the large-A^ theory. The running coupling in the Schrodinger 
functional (SF) scheme [98, 99] in SU(4) YM theory was studied in ref. [100], which found that, 
after converting to the MS scheme, the renormalized coupling nicely approaches the two-loop 
perturbative prediction as the momentum scale, at which it is defined, is increased. 

Motivated by the recent interest in walking technicolor theories [101, 102, 103], the authors 
of ref. [104] studied the mass anomalous dimension in generalizations of QCD based on gauge 
group SU(A'^) (with N = 2,3 and 4) with two dynamical flavors of fermions in the two-index sym- 
metric representation. Their results, displayed in the right panel of fig. 3, reveal clear similarities 
among the theories with a different number of colors. 

Lattice simulations of laige-N Yang-Mills theories at finite temperature T have been earned 
out in several works [105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116]: they revealed 
that all these theories have a finite-temperature deconfinement transition at a critical temperature Tc, 
which has a smooth dependence on (see the plot on the left panel of fig. 4, taken from ref. [116]). 
The transition is of first order for N >3, with a latent heat 0{N^) in the large-A'^ limit. 
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Figure 5: Left panel: Quenched mesonic spectrum at large A^, taken from ref. [126]. Right panel: Evidence 
for a rotor spectrum in large-A^ baryons of different spin, from quenched computations in QCD with N = 5 
and = 7 colors [127]. 



The equation of state has been investigated in various studies [109, 112, 114]: in the deconfined 
phase, the main equilibrium thermodynamic quantities, when normalized to the Stefan-Boltzmann 
Umit (which includes a A'^^ — 1 factor, related to the number of gluon degrees of freedom in the free 
theory) and plotted as a function of T /Tc, have almost no (residual) dependence on the number of 
colors A'^: see the plot on the right panel of fig. 4, taken from ref. [114], as an example. 

Similarly, good scaling properties with the number of colors hold for renormalized Polyakov 
loops in the deconfined phase, and their free energies for different representations are in excellent 
agreement with Casimir scaling [1 17], as recently shown by Mykkanen et al. in ref. [115]. 

The topological properties and the -dependence of large-A'^ gauge theories have been studied 
in various works [118, 119, 120, 121] — see also the review [122]. The topological susceptibility 
has a non- vanishing value in the 't Hooft limit, with small 0(1/ N^) corrections at finite A'^. 

Finally, there exist quenched computations of the large-A^ meson spectrum [123, 124, 125, 126] 
(see, e.g. , the plot on the left panel of fig. 5, taken from ref. [126]), and also of the baryon spectrum: 
in particular, the plot on the right panel of fig. 5, taken from ref. [127], shows clear evidence for a 
rotor spectrum (as expected from the arguments in refs. [128, 129]) in baryons of different spin in 
theories with N = 5 and N = 1 colors. 

2.2 Results in three spacetime dimensions 

Large -A'^ gauge theories in three spacetime dimensions (3D) share many qualitative features 
with their four-dimensional counterparts. In particular, the 't Hooft limit of these theories has been 
shown to be confining since many years [130], and, similarly to what happens in four dimensions, 
confining flux tubes can be described quite accurately in terms of Nambu-Goto strings, as shown in 
the left panel of fig. 6, taken from ref. [131]. Other recent studies of string effects in large-A^ gauge 
theories in three dimensions include works by Caselle et al. [132] and by Mykkanen [94]. 

In fact, string-like behavior has been observed in a variety of confining 3D models [133, 134, 
135, 136, 137, 138, 139, 140, 141, 142, 143] (see ref. [144] for a review): these studies confirm 
that the low-energy dynamics of confining flux tubes is consistent with the Nambu-Goto model 
at the lowest orders in an expansion around the long-string limit [145, 146, 147, 148, 149, 150], 
and are currently reaching levels of precision sufficient to reveal the deviations expected at high 
orders [151, 152, 153, 154, 155]. 
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Figure 6: Results in 3D YM theories: The left panel, from ref. [131], shows a comparison of the spectrum 
of confining flux tubes in SU(6) YM theory with the predictions of the Nambu-Goto effective string model. 
The right panel, taken from ref. [160], displays the trace of the energy-momentum tensor per gluon (in units 
of r^) as a function of T /Tc- The solid yellow line is the prediction from a 3D holographic model inspired 
by the improved holographic QCD model [161, 162, 163] (see also ref. [164]). 

For the YM spectrum, there is large numerical evidence that, like in four dimensions, glueball 
masses have smooth, finite large-A^ limits also in 3D [130, 156]. A very mild dependence on A'^ has 
also been observed for gluon and ghost propagators in Landau gauge [157]. 

Finally (and, again, similarly to the four-dimensional case), also the 3D equation of state 
appears to have only a trivial dependence on the number of colors [158, 159, 160], as shown, for 
example, by the plot of the trace of the energy-momentum tensor per gluon d.o.f. and in units of 
T^, as a function of T /Tc, on the right panel of fig. 6. 

2.3 Results in two spacetime dimensions 

In two spacetime dimensions (2D), several analytical (or semi-analytical) results have been 
known since the 1970's or early 1980's [165, 166, 167], but 2D theories continue to attract inter- 
est [168, 169, 170, 171] as useful QCD toy models. 

Recently, the eigenvalue density of Wilson loops in 2D has been investigated by Lohmayer, 
Neuberger and Wettig [172]; a related analysis in four dimensions has been carried out in ref. [173]. 

Finally, there exist studies of 2D large -A'^ QCD at finite chemical potential [174, 175, 176]. 

3. Concluding remarks 

Lattice studies of gauge theories in the \wge.-N limit are theoretically very appealing, numeri- 
cally tractable, and interesting for a very broad community. During the last fifteen years, numerical 
simulations in this field have given conclusive answers to various long-standing questions. How- 
ever, many other issues are still open, and waiting for your involvement. 

From my personal point of view, particularly promising research directions for further numer- 
ical studies at large A'^ include: 

• simulations with dynamical fermions, in various representations, which are potentially inter- 
esting also for physics beyond the Standard Model; 
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• further simulations at finite temperature and/or finite density, and comparisons with predic- 
tions from perturbative computations [177, 178, 179, 180, 181], from holography [182, 183], 
or from phenomenological models; 

• studies of the topological properties: although some of the open problems in this context are 
numerically challenging, the results of previous works are encouraging [118, 119, 120, 121, 
122], and there is steady algorithmic progress [184, 185]; 

• further investigation of \wge,-N equivalences and volume reduction: in the past few years, 
this subfield has seen a revival of interest, with impressive progress both on the theoretical 
and on the numerical side, and further works are well motivated. 
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